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^ Abstract 

iy\ We study the correction to the scale invariant power spectrum of a scalar field 

QQ on de Sitter space from small black holes that formed during a pre-inflationary 

£f) matter dominated era. The formation probability of such black holes is estimated 

P ° from primordial Gaussian density fluctuations. We determine the correction to 

£— s the spectrum of scalar cosmological perturbations from the Keldysh propagator 

(-***) of a massless scalar field on Schwarzschild-de Sitter space. Our results suggest 

— that the effect is strong enough to be tested - and possibly even ruled out - by 

observations. 
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1 Introduction 

Measurements of the cosmic microwave background (CMB) show that the Universe 
is isotropic on large scales to very good precision. Indeed, the CMB radiation is an 
almost perfect black body [1] with a temperature of T = 2.725 ±0.001 K 2 , and tiny 
temperature fluctuations superimposed with an amplitude of the order 10 _5 T |3]- The 
exception is the dipole, which is at the level of 10 _3 To, and can be explained by our 
motion with respect to the CMB rest frame. The assumption that measurements at 
any position in the Universe would lead to the same result implies that the Universe is 
also homogeneous. Observations of large scale structure support this assumption [4]. 
This Einstein's Cosmological Principle is a corner stone of modern cosmology and, 
hence, it is important to further experimentally test its validity. 

There has been much interest recently in the possibility that small violations of 
homogeneity and/or isotropy could give rise to the observed CMB anomalies. Some of 
the often quoted CMB anomalies are [6] : an anomalously small quadrupole and oc- 
tupole moments; a large deviation from the mean in some of the higher multipoles; the 
north-south asymmetry, the peculiar alignment of the quadrupole and octupole and 
their pointing in the direction of Virgo [7] ; the alignment of some of the higher mul- 
tipoles [5] ; the curious lack of power in the temperature angular correlation function 
on large angular scales [5], etc. Different authors disagree however in what consti- 
tutes significant deviation from homogeneity and isotropy. For example, Bennett et 
al. [5] tend to tune down the statistical significance of these anomalies, and argue that 
most of them can be attributed to priors. They also argue that, in the absence of a 
deep theoretical justification, which would make further tests possible, these anoma- 
lies will most likely remain curiosities. In addition to the CMB anomalies, there are 
also anomalies in the large scale structure of the Universe. For example, observational 
evidence was reported by Kashlinsky and collaborators [T2] for large scale (dark) flows 
of galactic clusters which cannot be explained by homogeneous, adiabatic, Gaussian, 
cosmological perturbations generated during inflation. Moreover, some authors [51110] 
offer an alternative to dark energy by considering the earth to be located near the 
center of a large void |10j , or by considering a randomly distributed collection of voids 
in the Universe (the Swiss cheese Universe) [TT] . 

The question we pose in this paper is whether some of these anomalies can be 
explained by placing a small black hole into an inflationary universe, r] Since an- 
swering this question rigorously is hard, here we make a first step in addressing it. 
In order to model cosmological perturbations we consider quantum fluctuations of a 
massless (or light) scalar field minimally coupled to gravity in Schwarzschild-de Sitter 
(SdS) space, and calculate the corresponding spectrum. Based on the knowledge of 
the Mukhanov- Sasaki gauge invariant potential, we then estimate the spectrum of 
scalar cosmological perturbations. We make the assumption that the gauge invariant 
treatment also applies to the inhomogeneous cosmology at hand, but warn the reader 
that our approach should be tested by a rigorous study of cosmological perturbations 
in inhomogeneous settings such as inflation endowed with a small black hole. In this 
paper we ignore tensor perturbations, since we expect that their amplitude will be, 
just as in homogeneous cosmologies, suppressed when compared to that of scalar per- 

3 Our study of small black holes can be quite easily extended to point-like particles such as magnetic 
monopoles or heavy particles whose mass is of the order the Planck mass. However, we expect that 
the effect of these particle-like objects will in general be much smaller than that of primordial black 
holes, giving thus a competitive edge to the study of black holes in inflation. 



turbations. Furthermore, the effect of vector modes is not taken into account because 
they are not dynamical for a homogeneous background and hence we expect them to 
play a subdominant role in the case of weak breaking of homogeneity that we consider. 

It is an important question how to relate our results for the primordial spectrum 
of scalar cosmological perturbations to the CMB observables. An interesting study in 
this direction is Ref . [T3] , where the authors investigate how different types of violation 
of homogeneity and isotropy would affect the temperature fluctuations in the CMB. 
Based on symmetry considerations, the authors consider in particular how a point-like 
defect (particle) , a line-like defect (cosmic string) or a plane-like defect (domain wall) 
would modify the observed temperature anisotropies. By symmetry a small black 
hole is closest to a point-like object, yet its event horizon makes it a more complex 
object to study, [j While the analysis in [T3j is useful to make a connection between 
cosmological perturbations generated in inflation and temperature fluctuations, it is 
not general enough to suit our needs. In particular, it cannot be used for a primordial 
black hole whose comoving distance from us is small in comparison to the wavelength 
of the perturbation considered. 

The main theoretical motivation for studying spectral inhomogeneities generated 
by a stationary black hole in inflation is that they yield results that can be tested 
against observations. This is so because the resulting spectrum can be viewed as a six 
parameter template. A good analogue are the gravitational wave templates provided 
by black hole binary systems. To illustrate more precisely what we mean, recall 
that homogeneous inflation produces a (power law) spectrum which, as a function 
of spatial momentum k, can be viewed as a two-parameter template, the parameters 
being the spectral amplitude (A^) and its slope (n s — 1), which have been by now 
tightly constrained by CMB measurements |14j . When viewed as a template, the 
SdS spectrum contains four additional parameters. These constitute the black hole 
position y with respect to us and its mass M during inflation r\ which we parametrize 
by fj, = (GMHq/2) 1 / 3 . Here Hq denotes the de Sitter Hubble parameter, and G is 
Newton's constant, ^j In the light of the upcoming CMB observatories, such as the 
Planck satellite, and ever increasing large scale galactic redshift surveys, it is clear 
that we will be able to test the inflationary black hole hypothesis. 

The paper is organized as follows. In section 2 we present the inhomogeneous 
Schwarzschild-de Sitter background metric and the corresponding scalar field equation 
of motion. We make a simple estimate of the number of black holes per Hubble 
volume at the beginning of inflation in section 3. Then, in section 4 we derive a 
formula which relates the amplitude of inflaton fluctuations on Schwarzschild-de Sitter 
space to the spectrum of scalar cosmological perturbations. The Schwinger-Keldysh 
propagator is derived in section 5 by expanding in the parameter fi. Section 6 deals 
with the application to cosmology. In particular, we obtain the power spectrum and 
illustrate its features in various plots in section 7. We close our paper in section 8 
with a discussion. Various technical details are relegated to four appendices. In 



4 A further complication is in the fact that a black hole could rotate and/or move with respect 
to the infiaton's rest frame. The latter could in principle be related to the claimed large scale dark 
flows [12] . 

5 Of course, the black hole has by now evaporated. 

6 If in addition the black hole is moving, three additional parameters are needed to specify its 
velocity; if it is rotating, three additional parameters are needed to specify its angular momentum; 
if it is charged, one more parameter is needed. We shall not study here observational consequences 
of these more general settings. Regarding the results presented in Ref. 1121 . it would be of particular 
interest to study the spectrum of a (slowly) moving black hole. 



this paper we work in units where H = 1 = c, but we keep Newton's constant G = 
6.674X1CT 11 m 3 kg- x s- 2 , the reduced Planck mass, M P = [SitG)- 1 / 2 = 2.4xl0 18 GeV 
and the Planck mass m v = G" 1 / 2 = 1.2 x 10 19 GeV. 



2 Inflaton field on Schwarzschild-de Sitter space 

2.1 Background metric and equation of motion 

A primordial black hole breaks the translational invariance of the background but does 
preserve rotational symmetry. The space-time metric of a black hole in an asymptot- 
ically homogeneous universe is the Schwarzschild-de Sitter (SdS) solution, giving rise 
to a line element which is usually written in static coordinates as 



ds 2 



-f(f)dt 2 



df 



,--' 



m 



f 2 dn 2 . 



(i) 



with /(f) = 1 — 2GM/f — Af 2 /3. In these coordinates the three symmetries of the 
SdS space are manifest: the time translation invariance and the two spatial rotations. 
Quantum fluctuations of a scalar field on the Schwarzschild background have been 
dealt with in Ref. [15] , where it was found that the radial mode functions of a massless 
scalar field can be expressed in terms of Heun's functions. But the presence of a 
cosmological horizon complicates the analysis and the SdS case has only been discussed 
for an extremal black hole [IB] • The reason for this is the difficult singularity structure 
of the d'Alembertian for (HI). For applications to cosmology another form of the metric 
is more useful. In Appendix A we show by explicit coordinate transformations that 
the metric takes the formQ 



ds z 



a z (r,)i -chf+[l + 



3 1\ 4/3 

/Try N 



1 - jj, 6 rf /r 6 

1 + /i 3 7y 3 /r 3 



dr A + r z dW 



(2) 



with [i = {GMHq/2) 1 / 3 and the scale factor a which is a simple function of conformal 
time 77, 0(77) = —l/(Hor)) (rj < 0). Notice that the metric ^ exhibits a black hole 
singularity at a finite radius, ro = — \xr\ (see also Eq. ( 108| ), such that in this metric 
r < r < 00 covers one half of SdS space. The Carter-Penrose diagram is plotted in 
Fig. fl] (see also [IB]) which also shows (schematically) how the interior of the black 
hole is covered by our coordinates. The Hubble rate Hq is related to the potential 
energy of the inflaton through the Fricdmann equation, Hq — V (<j>o) / (3Mp) . Here 
we assume that the inflaton potential energy V(4>o) is constant, such that it can be 
related to the effective cosmological constant as A = V((f>o)/Mp. 



The equation of motion for the massless inflaton field 
tion, 

D(j)(x) = 0, 



is the Klein-Gordon equa- 
(3) 



where the d'Alembertian □ acting on a scalar field is given by 

1 



□#a;) = ^VT^x) = -=d^V^gg^8 v <l>(x). 



(4) 



7 A similar form of the metric can be found in |17| . However, it has the disadvantage of being 
degenerate at the black hole horizon. 
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Figure 1: This plot shows the Carter- Penrose diagram of Schwarzschild-de Sitter space. It is an 
infinite sequence of black hole regions (Ilbh)i white hole regions (Ill w h), static regions (I and IV) 
and cosmological regions (II C and III C ). The coordinates we use reach the black hole singularity as 
well as the cosmological regions. We show schematically some r = const, lines (grey). In contrast to 
f = const, lines, they are timelike everywhere, cross the black hole horizon and eventually reach the 
singularity. Arrows indicate the flow in conformal time r\. Asymptotic future is given by r\ — > and 
corresponds to the boundary f = oo in the diagram. 



One easily finds from the determinant g of the metric tensor, 



g = a r sin 1 



pPrf 



(5) 



and hence 
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a 2 1- 
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where Vg 2 is the Laplacian on the 2-dimensional sphere, 



Vc2 = 



f a 



sin(0) 







1 d 2 



sm(9)d6 K '80 sin 2 (9) <90 2 



Because the CMB is highly isotropic, translation invariance in the early Universe 
can be only weakly broken [5, 6 . Hence /j< 1 and we can expand the metric in the 
parameter /x to first non-trivial order. The result is: 



ds z = a' 



drj 



' ^ 'l- 8 ££)dr> 



i+^v^^on+o^ 6 ) (7) 



3r 3 / \ 3r 3 

y/=g = a 4 r 2 sine? + O (ti 6 ) ; U = U ds + 5U + (^) . 

The differential operator D dS is the d'Alembertian on de Sitter space, 



° dS = ^ -^A+v 1 = A (-k 2 «+v 2 



a 
a 



(8) 



and 



with a' = da/dr] and V 2 = didi is the Cartesian Laplace operator in 3 dimensions. 

2.2 An estimate of the perturbation parameter \x 

The metric (|2j) contains the perturbation parameter \x = (GMHq/2) 1 ^ 3 as a constant. 



A more realistic point of view, however, is that a black hole of mass M decays due to 
an evaporation process and Hq is a time-dependent expansion rate of the Universe. 
Assuming a slow- roll inflationary scenario, i.e. small deceleration parameter 

e = -^«l, (10) 

we can neglect the time-dependence of Hq in the following discussion. Further com- 
ments on this approximation are made in section HI For the rate of change of M 
we make the following estimate. The evaporation time of a black hole is known to 
be PUHJ] 

T= 30720.^3 

9* 

where <?* ~ 10 2 — 10 3 is the number of rclativistic degrees of freedom at the energy 



scale ~ 1/(GM). The correction to the evaporation time (11) due to the Hubble 
horizon is negligible as long as \i <C 1 [20]. Assuming that the evaporation process 
takes longer than N = 60 e-foldings, we get 

'>-(l(S)f-. 

where we took GHq ~ 10 -12 and g* ~ 10 2 . In summary, for primordial black holes 
with a mass parameter 0.027 < /i « 1 we expect the SdS background to be a good 
realization of the inhomogeneous inflationary scenario. The constancy of \x is demon- 
strated in Fig. [2J 

3 Formation probability for black holes 

To assess the physical relevance of primordial black holes, we first estimate the prob- 
ability for their formation. For this we consider a pre-inflationary period dominated 
by heavy non-relativistic particles with a mass m < m p and Hubble rate iJ* < m p 
(here and henceforth an index * refers to quantities evaluated at the initial time £*). 
The presence of such particles is generally expected in physics at the GUT scale [21] . 
Moreover, such a scenario has recently been discussed in connection with CMB anoma- 
lies [2U H3] . On not too small scales the matter distribution in this pre-inflationary 
phase is well described by the local mass density p(x,t) — po{t)(l + S(x,t)). For 
black hole formation we are interested in the growth of density fluctuations inside of 
a bounded region. The number of particles in a comoving ball with physical radius 
R = R(t) and volume Vr at time t is denoted by N(R, t). The initial statistical fluctu- 
ations for the particle number are assumed to be Gaussian, 8N(R, t) ~ y / N(R,t) with 

G 
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Figure 2 : The dependence of /i on the number of e- foldings N is displayed here for different values 
Hq at the start of inflation. The red curve (second from below) corresponds to £t C rit ~ 0.027, the 
black hole evaporates after 60 c-foldings in this case. Below this critical value black holes evaporate 
before the end of inflation, black (lowest) curve. Above /^ cr it primordial black holes survive inflation. 
Their mass is seen to be constant to a very good approximation (blue, purple and brown curve). 



variance a(R,t) = (SN(R,t) 2 ) 
related to the local density, 



(N(R,t)). The mass density in the ball is linearly 



p(R,t) 



1 

Vr 



d 3 x p(x,t), 



(13) 



x\\<R 



which can also be written as p(R,t) = mN(R,t)/Vfi. Fluctuations in the particle 
number are thereby easily translated into density fluctuations. 

Black holes can potentially form from gravitational interaction of such fluctuations 
[15] . Due to the attractive nature of gravity, fluctuations can grow. To study their 
evolution, we consider the classical equation of motion for 5(x,t) = 5p(x,t)/po(t): r\ 

S + 2HS-c 2 J-j 6-4:TrGp {l + w){l + 3w)6 = 0, (14) 

with the speed of sound c 2 — (dpo/dpo)s — w, where po is the background fluid 
pressure. We take the equation of state parameter w — po/po to be constant. The 
Universe is assumed to be dominated by one species of particles in which case one 
can neglect entropy fluctuations, SS = 0, which would otherwise appear on the right 
hand side of Eq. ( 14 ) . When the gravitational term dominates in the above equation 
the density perturbation becomes unstable. The critical scale for the perturbation is 
given by the Jeans momentum 



k\ y/4nGpo(l + w)(l + 3w) 

a) j c s 



(15) 



which determines when thermal pressure is in balance with the gravitational force. 
The Jeans length A,/ = 2ir/(k/a)j reaches the Hubble scale if w « 1/3. For w « 
the Jeans length is very small compared to the Hubble radius and small black holes 
can form. We can solve Eq. (14) for k/a <C (k/a)j (super-Jeans scale) by making the 



Ansatz 8 oc t a and we find that 



8(k,t) = 8*(k) [ - 



■ (*)( r 



(16) 



3 c/. Bonometto (ed.), 'Modern Cosmology' (2002), p. 50. 



We shall neglect the second mode which is always decaying. In decelerating space- 
times (w > —1/3) the first mode is growing, whereas in accelerating space-times 
( — 1 < w < —1/3) fluctuations always decay due to the repulsive nature of gravity. 
The amplification actually increases with increasing w but, as mentioned before, we 
are only interested in the case w « since in this case \j/a <C Rh, which also has to 
be satisfied for the growing solution. The initial density perturbation is given by 

,,n = W) SN,(k) 1 

°* W ~ p m (k) ~ N*{k) ~ ,/NJk)' [U) 

Clearly, we have to look at large fluctuations away from the mean value (N) to find a 
significant probability of black hole formation. It is convenient to translate momentum 
space fluctuations Sp(k, t) to fluctuations in a ball of radius R by writing 

6N(R,t)_8p(R,t)_ r & MM)^), (18) 



(N(R,t)) Po (t) Jo ' po(i) 

with 5p(k,t) = (S p(k , t)) g ^ being the angle averaged momentum space mass density 

and 

w{kR) = l MkR)-kRcos(kR) 

■k kR 

is the well-known window function for spherically distributed matter in a ball of radius 



R. From relation (18 1 we conclude that SN(R,t)/(N(R,t)) grows precisely as the 



momentum space fluctuations (161, 



Next, note that 



SN(R,t) 

(N(R,t)) 


2(l + 3ra) 

5N*(R*) (H.A 3 < 1 +™' 


/ATI 


mlH 2 R 3 



(20) 



< 21 > 



which follows from the Friedmann equation, H 2 = (8ir/3)po(t)/m 2 and from 

3 (N(R,i))m 



Po(t) 



4tt R 3 



Now from Eq. (17) and (21) we can write, 



mlH 2 Rl 
<t.(R.) = (SN*{R,) 2 ) = (N*(R,)) = \ (22) 

Zm 



and hence 



4 ( 1 + 3 ") o „ 4 p6 - - - 4(1+3™) 



^ . ( ^ (fl ,t)) = <*w» 2 r^" 1 ^ = ^^ r^) 3(i+h . (23) 

V ^ V V '" (N*(R*)) \HJ 2mH 2 Rl\H ) V ; 

Due to the expansion of the Universe, the radius of the ball grows as R{t) — R* (a/ a* ) = 



i?*(i?*/-ff) 3(1+lu) in (23). In a decelerating universe, w > —1/3, the comoving radius 
grows slower than the Hubble radius, Rh = 1/H, such that, if i?* < 1/if* initially at 
t = t*, it will remain sub-Hubble at later times. This trend reverses during inflation, 
in which w < —1/3. 



A black hole with Schwarzschild radius Rs forms if, due to statistical fluctuations, 
the number of particles in Vr s becomes sufficiently large, N(Rs,t) > Rsrrii/(2m), 



(m 2 p = 1/G). Using ([21] and writing N(R s ,t) = (N(R s ,t)) +SN(R s ,t) one obtains 
the condition 



mzRq raiRa 

5N > SN cr = -^ - (N(R s ,t)) = -*-?-[l - (HRs) 2 } 



2m 2m 

Thus, using that the fluctuations SN are Gaussian distributed Q 



(24) 



P(SN) = 



1 



=cxp 



5N< 



J2-KG " V 2a 
with a = (SN 2 ), the probability that a black hole forms is found to beFj 

/ SN 2 

I" 00 1 CX P _ ?rr(R 

P(SN(R s ,t) > 6N CI (R s ,t)) = / d(SN)P(SN) « —= V ^ 

j(5AT cr ^v^ oJy cr /V2a 



(25) 



where, making use of Eqs. (23 1 and (24) 



<5iV r - 



y/2a(R,t) 2H*y/mR* 



[l-{RHf] 



H 



2H\ 



[l-(^) 2 ](f) 3<1+ "» 1.(26) 



The inequality in (261 is needed to correctly evaluate the integral (25). Notice also 



that, when that inequality is met, the probability for black hole formation is (exponen- 
tially) suppressed. Clearly, the inequality is broken for super-Hubble scales, for which 
RH > f . But at super-Hubble scales we expect suppressed statistical fluctuations, and 
hence do not trust our analysis anyway (see the comment further below) . Remarkably, 
up to the f — (RH) 2 term, Eq. (26) is time independent, which also means that the 



probability for black hole formation ( 25 1 in a ball of constant comoving radius R is 



time independent. Hence, the growth of perturbations (20) precisely compensates the 



decay in fluctuations (17 1. This might be telling us something deep about gravity. 



However, we do not have a simple explanation for this fact. For later purposes it is 



useful to rewrite our result for the probability of black hole formation (25 26) as 



P(n,m,H) = 



/' ; 



3/2 



lAmH 



1 - 16/u 6 






exp 



ml(\ - 16//) 



6\2 



H 



lGmHfj, 3 

(27) 
where we made use of the mass parameter /x = (GMH/2) 1 / 3 . Obviously, /i < 4 ^ 
must be satisfied in order for a black hole to be sub-Hubble. The question is then 



how to convert the probability ( 25 ) into the number of black holes formed before 



inflation and during inflation. The analysis presented above is meant to provide a 
rough estimate of the number of sub-Hubble black holes formed before inflation, and 
neither takes a proper account of causality, nor of nonlinear dynamics of over-densities. 



9 Recall that, from the central limit theorem, the Gaussian distribution is the large N limit of the 
Poisson distribution. 

10 In the limit SN 2 2> a the statistical fluctuations could also obey a power law behavior, 
■fcrit(^-N') a (u/SN 2 ) x , like in the theory of critical phenomena. We are not going to consider 
this possibility here in any detail. Note, however, that for this type of statistical fluctuations more 
black holes will form. 




Figure 3: In this figure we display the formation probability of black holes per Hubble volume on a 
logarithmic scale as a function of the mass parameter [i = (GMH/2) 1 ' 3 for particle mass m = 0.3m p 
and initial Hubble parameter H, = 0.3m p (solid curve, red); m = 0.3m p and H* = 0.1m p (dashed, 
blue) and m = 0.01m p and H* = 0.3m p (dot-dashed, black). The left vertical dashed line indicates 
the critical black hole mass /x cr it = 0.03 and the right one marks fi = 4 — 1 ' 3 . Black holes with mass 
parameter close to /^ cr it are seen to be Gaussian suppressed, whereas in the intermediate range their 
number can be Oil) per Hubble volume at the beginning of inflation. 



Staying within this type of reasoning, we propose to interpret (27 1 as an estimate for 



the probability that a black hole formed by the beginning of inflation in a comoving 
volume Vr — (47r/3)i? 3 . The expected number of black holes per Hubble volume 
(4tt/3)R^ 1 at the beginning of inflation is then 



(N Ba {n,m,H)) 



P(l*,m,H) 

(RsH)z 



P(»,m,H) 
64/x 9 



(28) 



If the Hubble volume today corresponds to X Hubble volumes at the beginning of 
inflation, then there will be about X(7Vbh(m, fn, H)) pre-inflationary black holes within 
our past lightcone. A more detailed discussion on how to relate the number of black 



holes (28) to the number of (pre-)inflationary black holes potentially observable today 
is given in section [8J In Fig. [3] we show how the expected number of black holes ( 28 ) 
depends on the black hole mass parameter /i for different values of particle mass m 
and initial Hubble rate H*. We emphasize, however, that in our analysis we make 
the assumption that the statistical fluctuations are normally distributed on all scales. 
This might not be true, as has been argued for example in [23], where fluctuations 
of other types are considered on super-Hubble scales. Causality can limit the size of 
these statistical fluctuations. It may be more realistic to assume that on super-Hubble 
scales surface fluctuations are dominant, SN oc y/S ~ iV 1 / 3 and a w (N) 2 / 3 , thus, 
suppressing the formation of black holes that are initially super-Hubble. The actual 
probability for black hole formation might therefore be smaller in the region fi ~ 4~ 1 ' 3 
than it is shown in Fig. [3j Moreover, the formation probability depends strongly on 
the mass m of the heavy particles, and yet we do not know much about it. Based on 
our current understanding of particle physics and gravity, it is reasonable to assume 
that m is limited from above by the Planck mass m p . If to ^C m p , however, the 
particles would start behaving relativistically, which would increase the Jeans length 
and further suppress, or even prevent, the formation of black holes. 
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4 From scalar fluctuations to scalar cosmological per- 
turbations 

4.1 Homogeneous background 

Before discussing the effect of inhomogeneities on the scalar spectrum, we review the 
treatment of cosmological perturbations on homogeneous backgrounds, such as the 
conformally flat background metric, g b = a 2 (7/) 77^, where rj^ v — diag(— 1, 1, 1, 1) and 
a is the scale factor of the Universe. The physical situation we have in mind is a slow- 
roll inflationary model driven by a homogeneous inflaton field <po(t), where the Hubble 
parameter H (t) = a/a (a = da/dt = da/(adr])) is a slowly varying function of time, 
such that the de Sitter limit is obtained when H — > H = const. Scalar perturbations 
in this model are induced by the quantum fluctuations of the inflaton, while tensor 
perturbations are induced by the quantum fluctuations of the graviton. In linearized 
perturbation theory the two decouple. For a recent treatment we refer to |25j . while 
standard reviews are Refs. [261127], 

It is convenient to decompose the inflaton $ and the metric tensor g^ into the 
background fields 4>o,gtv an d the fluctuations (/>, /i M „ as follows, 

#(s) = Mt) + 4>{x) , g^{x) = gl„(t) + a 2 h^{x) . 

The background fields a and <po are classical field configurations, whereas hnu and <fr are 
dynamical quantum fields. A detailed study (see e.g. |25j ) shows that there are only 
three physical degrees of freedom, two from the graviton and one from the scalar field. 
These can be expressed in terms of the gauge invariant Mukhanov-Sasaki variable 
v = —((f>o/H)lZ (4>q = d(f>o/dt), and the gauge invariant tensor hf? , where by gauge 
invariance we mean the invariance under linear coordinate shifts x^ — > x^ + ^{x). 
Since physical observables are independent of the choice of gauge, it is convenient to 
work with gauge invariant fields. In homogeneous cosmology the Sasaki-Mukhanov 
field (curvature perturbation) is of the form, 

H = i>-^4>, (29) 

where ip is the scalar gravitational potential defined by the scalar-vector-tensor de- 
composition of hij : 

ha = 2iP6 %3 - IdidjE + 2d {l F }) + h^ (30) 

and (j> is the inflaton fluctuation. The potential ip is a gauge variant measure for local 
spatial volume fluctuations. When working with gauge invariant variables, such as 



1Z in (29 1, we are guaranteed to get observable CMB temperature fluctuations, since 
it is the gradient of 1Z that sources photon number fluctuations through the photon 
Boltzmann equation. One can also get a physically sensible answer when one fixes 
a gauge, provided one makes the correct link to the late time gravitational potential 
which enters the fluid equations. 

For example, in the comoving gauge, in which <f> = = E = Fi = 0, it is the spatial 



gravitational potential ip that determines the Sasaki-Mukhanov field, TZ = ip (29). On 



the other hand, in the zero curvature gauge, in which ip = = E = Fi = 0, it is the 



inflaton fluctuation <f> that determines TZ through ( 29 1 



TT 

TZ = — —(f> (zero curvature gauge). (31) 

00 



If 



This relation can be used to estimate the late time potential from the inflaton fluctu- 
ations 4>(x). The spectrum of scalar cosmological perturbations V-r. in zero curvature 
gauge, (V-r. is conserved on super-Hubble scales) is related to the spectrum of scalar 
field fluctuations V^ as, 

H 2 „ 1 

2eM'i 



Vk = ^V^ = —^V^ (32) 



'u 



L p 



where M P = I/V&kG = 2.43 x 10 18 GeV = 4.34/Ltg denotes the reduced Planck mass. 
In order to get the latter identity, we used the second Friedmann equation, —4:irG(j> 2 — 
H, and the definition for the slow-roll parameter, e = —H/H 2 . 

Of course, there are two remaining degrees of freedom in the graviton which have 
not been taken into account. However, the graviton spectrum V grav is known to be 
suppressed in the slow-roll inflation [28] HH [25] 



*W=^, -^ = 16e«l, (33) 

such that, to first approximation, we can neglect the graviton contribution to the 
spectrum of cosmological perturbations. 

Finally, the field fluctuations can be translated to the temperature-temperature 
correlation function as [T3"j . 

ST(h 1 )6T{h 2 )\ _ fd 3 kd 3 k' ^(2l + l)(2l> + l) r l+l , 



E 'fZv' W'** &> 



n I j w jjf (4tt 

x P t (k ■ hi/k)Pi>(k' ■ n 2 /fc')^(fc)6 i (A?)^(fc')e/'(fc') , 

where Pi{x) denotes a Legendre polynomial, and Qi(k) and 0j/(fc') denote the ap- 
propriate transfer functions, which relate IZ(k) to 5T(k), and which are obtained by 
solving the Boltzmann equation for the photon fluid. 

4.2 A small black hole in a de Sitter Universe 

Keeping in mind the procedure of the previous section for the derivation of cosmo- 
logical perturbations from scalar fluctuations, we have to pay special attention to use 
gauge invariant fields also in the inhomogeneous case. We split the metric into the 
Schwarzschild-de Sitter background (17]) and a perturbation h^ v . The spatial part of 
the perturbation hij can be decomposed as 

ha = 2V>7ij - 2 V (i V j) E + 2V {t F 3) + h T J , (35) 

where jij is the induced metric on the spatial slices for the metric (M) and V is the 
covariant derivative compatible with jij. The vector £ M that generates the gauge 

transformation x M — > x^ + £ M (x) can be written as £ M = (£°,£ I ), with £* = £p + V £ 
and Vj£r = 0. The metric perturbation transforms as 

hij -> h'ij = hij - 2V (i ^) (36) 

= 2ihii - 2V (i V i) (E + O + 2V (i (f>) - Q + h% T + 2r?^o , 
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with the Christoffel symbols of the metric (frl) 



^ = % - ^*r« (37) 



a " J 2a 



and 7ij = 7?? + 5~/ tj + 0{p 6 ) with 



10 \ ,0,3 / -2 



7^ = a 2 ( r 2 I , 5 7ij = ^f^- (Or 2 | . (38) 

r 2 sin 2 9 I V r 2 sin 2 



4a n A ri 

3r 3 

Moreover, we observe that 5jij is traceless to the relevant order in /i, showing that £0 
generates also transformations of Fj and hf T . The transformation properties of the 
two scalars ip and E are 

■<Jj^iP+—£o+0(p e ), E^E + Z + Obi 6 ). (39) 

a 

Together with the gauge transformation of the scalar field, <fi — >• (f>+ (0o/° 2 )Co, WG hud 
that 

V-^^-^ + ^M 6 )- (40) 



We conclude that the Sasaki-Mukhanov field 1Z in Eq. ( 29 1 is gauge invariant (to order 
/U 3 ) also when the inhomogeneity caused by a small black hole is taken into account. 

The fact that £0 generates gauge transformations of a vector and the tensor sug- 
gests a mixing of the scalar, vector and tensor sectors in inhomogcneous cosmology. 
Therefore, for a completely rigorous treatment one should look at the quadratic action 
on the SdS background and take into account the couplings of the modes. 

Furthermore, we make use of the slow-roll paradigm. This means that, even though 
strictly speaking our results will be derived in SdS space, we shall assume that they 
hold in quasi de Sitter space endowed with a small (decaying) black hole, provided 
one exacts the replacements: H —> Hit) and M — > M(t). This is justified when both 
the Hubble parameter and the black hole mass change adiabatically in time, in the 
sense that H < H% (e < 1) and M < MH . 

5 The propagators 

To study scattering experiments, one typically calculates the S-matrix elements. In 
cosmology, on the other hand, one is primarily interested in expectation values of 
operators with respect to some definite vacuum state. For this purpose the in-in, or 
Schwinger-Keldysh [211 [30], formalism [3"Tl I521 12"5] is suitable, in which time evolution 
of an operator is described in terms of the perturbation theory based on the Keldysh 
propagator and the in-in vertices. Since here we are primarily interested in the 
spectrum of a scalar field on SdS space, which can be obtained from any equal time 
two-point correlator, for our purpose it suffices to calculate the corresponding Keldysh 
propagator. 

The Keldysh propagator is a 2 x 2-matrix of the form, 

iG{ X ;x')={ iG r ++ f ]X 'l iG r + 'f ]X '\), (41) 
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whose components are the Wightman functions iG^ , iG |_ and (anti-)time ordered 

Feynman propagators iG ++ , iG , defined as, 

iG F (x;x') = iG ++ (x;x') = (Q\T(j)(x)4>{x')\n) (Feynman) (42) 

iG+-(x;x r ) = (n\c/)(x')<j>(x)\n) (Wightman) 

iG-+(x;x') = (n\<j>(x)<f>(x')\Q.) (Wightman) 

iG—{x;x') = (n\T(j>(x)4>(x')\Q,) (anti - Feynman) , 

where |fi) is a suitably chosen vacuum state. The time ordering is defined as 

TcJ)(x)<j)(x') = 9(x — x' )cf>(x)<f){x') + 9(x' — x )4>{x)(j)(x) (43) 

T4>{x)4>{x r ) = 6(x - x' )cf>(x')<j>(x) + 9{x' - x )<j>(x)cf>(x') , 

i.e. later times are to the left for T and early times are to the left for T. The propagator 
iG satisfies the equation 



\/-g(x) (□. - <) iG{x- x') = inH\x - x') , (44) 

where <r 3 is the Pauli matrix 

° 3 = {l ->) <«> 

and m 2 s = d 2 V((f>o)/d<j>Q is the effective mass-squared of the field, which in the fol- 
lowing we neglect. In slow-roll inflation m a s can be expressed in terms of the second 
slow-roll parameter T)y — MpV"/V, i.e. m 2 s — r]yV((f>o)/Afp, such that setting 
"T-eff — >■ is equivalent to r\y — > 0. 

5.1 The de Sitter case 

In the de Sitter case we can solve the equation of motion for the massless scalar 
field (pi [4]) explicitly. Taking advantage of spatial homogeneity of de Sitter space, the 
following mode decomposition of the free field is convenient, 

<Kv, *) = / (S3 [«* (V)^ k - S b(k ) + u* k fa)e-**-*&t (k )] , (46) 

where b(k) and b'(k) are the annihilation and creation operators, defined by b(k)\il) = 
and by the commutation relation, [b(k), b*(k')] = (27r) 3 5 3 (/c — k'). The mode 



functions u^ (77) in ( 46 I satisfy the equation 

2 



d 2 n + k 2 - — 1 {au k {rf)\ = . (47) 

We obtained this result by making use of ^ and noting that a(rj) = —l/(H r]) implies 
a" /a = 2/rj 2 . Imposing the boundary condition that the mode functions behave like 
in the conformal vacuum in the asymptotic past yields 

UkM = --^=(l-^\ e- ifc " . (48) 

U ' a^2k\ k V J v ; 



This equation, together with the condition b(k)\fi) — (for all k), defines the Bunch- 
Davies (BD) vacuum |17). The fact that in the ultraviolet (k/a 3> H ) the BD vacuum 
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minimizes the energy in the field fluctuations has led to the belief that this vacuum 
represents a sensible physical choice [33] for the inflationary vacuum. This is not 
so because in the infrared (where k/a <C Ho) the BD vacuum yields infinite energy. 
Namely, in the IR adiabaticity of the state is broken because the field couples strongly 
to the expanding background, leading to abundant particle generation, having as a 
consequence strongly enhanced infrared correlations. While this particle creation is 
very welcome in cosmology, since the amplified vacuum fluctuations provide a beautiful 
explanation for the Universe's structure formation, one has to take proper care to 
regulate the IR. One way of doing that is to replace the BD vacuum by a more general 



state, characterized by the following generalization of the mode functions ( 48 ) 



u k (r))^u k (r)) = a(k)u k (n) + /3(k)<(v); \a(k)\ 2 - \p(k)\ 2 = 1 . (49) 

By suitably choosing /3(fc), one can then make the infrared part of the vacuum state 
finite 34 . A concrete working realization of this proposal has been investigated in 
Refs. [35, 36 . Alternatively, one can remove the infrared problems by placing the 
Universe in a large comoving box of size L. This leads to a discretized reciprocal 
(momentum) space k = Tifco [ft — (711,712,713), with rii integers], with the comoving 
lattice size k = 2tt / L. Since fco corresponds to the minimum allowed momentum, this 
cures the infrared problem simply by disallowing the deeply infrared modes. In the 
limit when L — > 00, the lattice constant fco — > 0, and the sum over the momenta can 
be replaced with increasing accuracy by an integral, which has fco as the IR cut-off, 
thus regulating the infrared. From the physical point of view, it is natural to associate 
this cut-off with the scale of the Hubble horizon at the beginning of inflation, thereby 
eliminating modes that stretch beyond the Hubble radius. One way of implementing 
this, is to take the spatial topology of the universe to be compact, e.g. a torus, as 
discussed in [57) . 

To see how the regularization procedure works in practice, we shall now calculate 
the regulated Feynman propagator in de Sitter space. In order to do that, we need to 



relate the direct space propagator to its mode functions ( 48 ) . Because de Sitter space 



is spatially homogeneous, it is convenient to write the components iA a i>(x; x') of the 



Keldysh propagator (41 1 for de Sitter space in terms of its Fourier space counterparts, 



1 rfik - / 

iA ab (x; x') = / ^e <k -<*-* ) iA afc (fc, 77, 17') . (50) 



Making use of Eqs. mm and (48), one finds for the momentum space propagators, 

iA + _(k,r),ri') = u* k {rj)u k (r)') 

*A_ + (fc,77,r/) = Uk(v)u* k (r]') = (iA + _(fc,?7,r/))* 

iA++ (fc, 77, 77') = 6(77 - r/)iA.+ (fc, 77, 77') + @(r)' - 77)iA+_ (fc, 77, 77') (51) 

iA__ (fc, 77, 77') = 6(77 - 7/)iA + _ (fc, n,rf) + 6(77' - ?/)iA_+ (fc, 77, 77') . 

The corresponding spectrum Vi , defined by 

m(s,v)4>(s , ,vm = jfrf ^ ^^kW^f 1 ' (52) 

is obtained straightforwardly from the equal time limit (77' — > 77) of the propagator, 
Vf (M) = |^A + _(fc,77,77) = |j| (1 + fcV) ■ (53) 
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It is scale invariant at future infinity, 77 — ^ 0. 



Based on (51) one can calculate the position space de Sitter propagator by per- 



forming the momentum integral (50) over k > k$. The resulting Feynman propagator 



is 



iAp(x; x') = 



Hi 

Ait 2 



Ax 2 



log (k 2 Ax 2 ) + 2( 7E - 1) + O(k ) 



(54) 



where Ax 2 = — (|ry — rf\ — it) 2 + \\x — x! || 2 is the conformal space distance function. 
Two comments are in order. Firstly, apart from the standard Hadamard contribution 
oc 1/Ax 2 , which is singular on the lightcone (on-shell) and quickly decays off-shell, due 



to rapid particle production in de Sitter space, the de Sitter propagator ( 54 ) acquires 



a logarithmic term which contributes both within the past and future light cones. 
Secondly, the logarithm grows without a limit as k§ — > 0. This is a manifestation of 
the IR singularity of the Bunch-Davies vacuum. We will see below that a black hole 
in de Sitter space 'sees' this logarithmic singularity in the corrected SdS propagator 
as a logarithmic singularity in the mixed space propagator and hence also in the 
(mixed space) SdS spectrum. This dependence on the IR regulator poses a unique 
opportunity to investigate the black hole contribution to the spectrum dependent on 
the IR regularization. In this paper we choose the comoving box regulator primarily 
because of its simple implementation, but one would certainly benefit from studying 
other IR regularization schemes. 

5.2 The Schwarzschild-de Sitter case 

For the case when a primordial black hole is present in de Sitter space we shall derive 
only the first order correction in fi to the Schwinger-Keldysh propagator. For this we 
write 



□ = □ 



dS 



SO: 



/ -0 = V-Qds + $V-g ■ 



iG = iA + iSG , 



(55) 



into (44) we find that the correction to iSG satisfies: 



with iA being the propagator on de Sitter space (|50f51), (54). By plugging this 

(56) 



V^g^M (Df - m c 2 ff ) tSGix; x') = 

- \/-9ds{x)Sn x iA(x; x) 



<V-s( x ) 



a 3 8 i (x - x') 



V-9ds(x) 
Note that Sy^g is only C(/x 6 ), and we will neglect it from now on. It follows that 

i5G(x] x')=i J d A x"yJ-g ds {x")iA{x] x")a 3 8D x „iA{x"; x') . (57) 



This solution of ( 56 ) is given only up to a homogeneous solution of the d' Alembertian 



operator in (56). The unique propagator in (57) is obtained upon specifying the 



boundary conditions for the mode functions, or equivalently, for the vacuum state. 
Here the unperturbed vacuum state is chosen to be the (pure) Bunch-Davies vacuum 
of de Sitter space, whereby the deep infrared modes are removed by placing the Uni- 
verse in a comoving box, as explained in section |5.1| But we are still free to add a 
homogeneous solution to the Feynman propagator (resulting in a mixed state) , which 
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has to be added to take Hawking radiation into account. The light black holes that 
we consider do indeed emit Hawking radiation but a simple estimate P_H suggests that 
this does not change the spectrum of scalar fluctuations at the leading order in the 
perturbation parameter /i. Namely, comparing the emission rate per Hubble time tjj 
and Hubble surface area An for inflaton fluctuations and Hawking radiation, we find 
that 



E 



A u t 



H L H 



m 



E 



A„t 



H L H 



MHl 



Rs_ 
R H 



H?H*, (58) 

Haw. rad. 

where the factor (Rs/Rh) 4 accounts for the redshift of Hawking radiation from the 
time of its creation, when the typical wavelength is of the order of the Schwarzschild 
radius, to the time when the amplitude freezes out, when the wavelength is of the 
order of the Hubble radius. This estimate holds in the case when the scalar field 
is light (or massless), m <C Hq. Note that there is no enhancement in the second 



equation of (58 1 by the total number of degrees of freedom g* (cf. Eq. ( 11 )) because 
only non-conformally coupled matter fields (scalars and tensors) are relevant for the 
radiation from the black hole as their amplitude freezes out after exiting the Hubble 
radius. From |58|) we find that the ratio of the emission rates is suppressed as ~ /i 6 



since M ~ M%/M 2 oc ^ b and (Rs/Rh) 4 



I' 



The emission of very massive 



particles, 



< GM, is exponentially suppressed, making them irrelevant for the 



above estimate. In the intermediate regime, GM < to < Hq , the ratio of the 
emission rates is ~ ^(Hq/wl). This dimensional analysis indicates that Hawking 
radiation does not contribute at the leading order which is \x or /i 3 , depending on the 



scale, as can be seen in (74) and (71) arising from the homogeneous contribution to 
the propagator. 

Since i5G(x; x') — > for M — > 0, our modified vacuum state reduces to the Bunch- 
Davies vacuum of de Sitter space in this limit, as it should. When Eq. (57) is written 
in its component form we get, 

i5G ab (x;x') = i J2 c J d i x"^~g dS (x")iA ac (x-x")5n x „iA cb (x"-x l ) 1 (59) 



with a, b = +, — . Writing the propagator iA in momentum space ( 50 51 ) this becomes 
i6G ab (x]x') 

4i[i 3 [ d 3 kd 3 k' f" 

coo ^ r ii r 1 r 27r I 1 1' . -"V-' 



(60) 



W E c l A ac (k,r 1 ,7 1 ") l A cb (k / ,7f,r 1 / ))e^ s - S '^ 



dcos9" / d<t>" -fc' 2 + 3 



i (jj ' x ) \ i(k'-k)-x" 



..112 



It turns out that it is easiest to evaluate a (double) momentum space version of iSG ab - 
To do that, we first introduce the momenta associated with the positions x and x' , 



i5G ab (f u p 2 ,r],ri') = / d 3 x d 3 x' iSG ab (x,x' ,??,r/)e ifl ' 



(61) 



11 From Table 1 in 1151 one can in principle obtain an explicit expression for the asymptotic form of 
the (renormalizcd) propagator for Hawking radiation in the Unruh vacuum, which is the physically 
relevant state in this case. However, a rigorous treatment of the effect of Hawking radiation requires 
the knowledge of the asymptotic solutions to the confluent Heun's equation which determine the 
radial mode functions but these asymptotic solutions are, to our knowledge, not known. A detailed 
analysis which could support our conjecture in Eq. 158|| is beyond the scope of this work. 
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and, next, relative and average coordinates in position and momentum space, r 
x — x', y — (x + x')/2 and p = pi + P2, q = (pi — P2)/2. This yields 



iSG ab (p,q,r,, V ') = J d 3 yd 3 ri5G ab (x 7 x' ,77, t/K^V^ 

,3 



(62) 



= - 3H 2 fc ,2 (3cOS 2 6 - 1) ^ c ( J ac,cb{^P^^V,v')- Jac,cb(Vo;P,q,V,v')) ) 



c=+,- 



where here and in what follows k = q+ \p, k' = q — - 
<^( — PtQ^ \p)- Moreover, we defined 



i" ' : ' - * - \p and 6 = <(k' - k,k') 



/ ,, ,s f cosjupr]") _ smjfiprj") 



JabMv";P,q,V,v')= dr,''r,''i^ ah {k,r,,r ) ")i^ cd {k',r,",r ] ') ™^ 
J \ [i l p A r\ 



Y~i' 1 £ ll^ Yi nri'*& 



(63) 

The details of the derivation can be found in Appendix B. From Eq. (60) and rela- 
tions (111 I for the step functions, we find that the corrections to the Feynman and 
anti-Feynman propagators obey the standard time-ordering relations 

iSG ++ (p, q, 77, 7/) = 9(77 - r]')iSG- + (p, q, 77, 77') + 9(?/ - r])i5G + -(p, q, 77, 77') (64) 
iSG — (p,q,r],r)') = 9(77 - i]')i8G+-(p, q,r), 77') + 9(77' - rj)i6G-+(p, q, 77,77') . (65) 

In addition, we have 



i8G- + (p,q,r),r)') = (i5G + -(p,q,ri,r)')) 



(66) 



Therefore, in order to fully reconstruct the black-hole-corrected Keldysh propaga- 
tor,we only have to determine iSG^ , for which we need to know only J_| _| and 



J^^M^ ^f^ (i + ±) (1- JL) e -'"-'V W) 



4kk'(np) 4 \ krj 

(^) 2 7( M p) 2 (fc~fcQ77"~<# 2W 

kk' (k- k') 2 - (jip) 2 ((k - k') 2 - (up) 2 ) 2 



I 

k'rf 

\4 



W (M 3 77"+ t( ^- fc ' ) +7(^)(fc-fcQ 2 l ( W f(k-k>) 

kk'r," (k - k') 2 - (pp) 2 ((k-k 1 ) 2 - (up) 2 ) 2 



cos([ipr/") 

sm(npr]") 



H^nn'e i< - kv+k ' n '^ I 1 

J+--+(V ,P,Q,r),V) = rr^r; — I I 



4kk'()ip) 4 

(w? *(w) 2 (k + kiWi + ^ 
kk' 



1 



krjj \ k'rj' 

2(M 4 



(k + k 1 ) 2 - (up) 2 ((k+k') 2 - (fi P ) 2 y 



-i(k-\-k )r] 



cos(fj,pr] ,f ) 



(68) 



HP , (wOW - lM l^' ] + i(w>)(fc + fc') 27( w ) 3 (fc + fc') 



fcfc'77'' 



(A; + fc') 2 - (mp) 2 



((fc + fc') 2 - (/ip) 2 ) 2 



sm.(p,prf') 



18 



From Eq. ( 62 ) we get 



i5G+- (p, q, n, n') 



167TI/1 3 



k' 2 (3cos 2 8~ l)[j;_ i _ + ( I7 ) - J+-,+-(v) 



(69) 



+ J+-.+-W) - J+-.-+W) + J+-,-+(vo) - Jl-,-+Cno)] 



where here, for brevity, we wrote J a b,cd(v") — Jab,cd{ , n"]p,q,rj,ri l ). In order to get 
the complete SdS propagator, we still need to add the de Sitter propagator to the 



correction (69 1, which in the double Fourier space reads, 



HgiW)' 



iA + _(f,q,r,,r,') = ^-(l + -) (l - —V^-^) 3 ^) • (70) 
2a V an/ V 077'/ 



1 + — 
qrj 



1 

qrf 



Together with Eqs. ( 64 67 1 and (70), relation (69) completely determines the desired 



SdS propagator in the limit of a small black hole mass, and it constitutes our main 
result. An interesting feature of the propagator ( [69] ) is that, due to causality, it does 
not contain any information from future infinity (77" = 0). Notice that, although 
/i <C 1, we have not expanded the factors sin(/ip?/') and cos(/ip?/') in powers of [i 
in Eqs. (67-68) since one might still want to consider momenta p with —fiprj" > 1. 
Consequently, even though our original expansion parameter was /u 3 , the propagator 
correction is formally suppressed only as fj, 1 ex M 1 / 3 , thus, as a fractional power of 



the black hole mass. Finally, it is worth noting that the pole of J+ 



at k + k' 



fip 



in (68) is not physically realized for any /j, < 1, because fc + fc' 



||g+p/2|| + ||g-p/2||> 



6 The power spectrum 

6.1 Double momentum space representation 



In section 5. 1 wc derived the spectrum ( 53 1 for a massless scalar field on de Sitter space. 



The inhomogeneous case with a primordial black hole is far less trivial to deal with, 
mainly because out of the 10 symmetries (Killing vectors) of de Sitter space only three 
symmetries remain in Schwarzschild-de Sitter space. (Recall that the homogeneous 
cosmology of slow-roll inflation, radiation and matter era has six symmetries.) 

In principle, one has to rederive the gauge invariant combinations of the fields, 



such as the Sasaki-Mukhanov field (29) TZ, for the inhomogeneous background and 



determine their power spectra. But we have seen in Eq. ( 40 1 that TZ is approximately 



gauge invariant for weak breaking of translational symmetry. Similarly, the graviton 
contribution can be neglected, since its power spectrum is expected to be equally 



suppressed, Eq. (33) 



With these assumptions the correction to the spectrum from the black hole can be 
determined from i5G + - by taking the equal time limit n' -> 77 of the propagator ( 69 ) , 



i5G+- {p, q, 77, 77) 



327T/i 3 



fc'^(3cos^6»-l) 



(71) 



[Im J + _ t _ + (r) ;p, q, 77, 77) - Im J + _,_ + (?7; p, q, 77, 77)] 



The 77 dependence of the corrected power spectrum SV ! p(p, q, 77) is displayed in Fig. El 
Note that, because of Eq. (32), the relative correction to the spectrum of inflaton 
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Ho 3 ^,. (17) 




1 Ho 



HJ5P4, (17) 

P/ S (V) 




II l'„ 



Figure 4: The dependence on r) of the correction to the power spectrum in (double) momentum 

space, SV l / > (p,q,ri) = (q 3 /(2n 2 ))iSG^ (p,q,r],r]), rescaled by the de Sitter spectrum "P<p(<f, ??) = 

(Hq/(4tt 2 ))(1 + q 2 r] 2 ), is presented here for fixed momenta p, q and different values of fi (from 
bottom to top: fi = 0.1 (blue curve), fi = 0.09 (pink), /j = 0.08 (yellow) and fi = 0.05 (green)). The 
correction is zero at the initial hypersurface, r\ = —Hq . For small r\ the spectrum is well described 
by an expansion to order rj 3 . The linear order vanishes. It approaches a non-zero value in the limit 
T) — > 0. Left panel: p/Ho = q/Ho = 1 and <(p,qj = 0, right panel: p/Ho = 10, q/Ho = 1/10 and 
<(p,q) = n/2. Note that in this and in the subsequent plots the relative correction to V^ is shown 
but this coincides with the relative correction to Viz because of Eq. l|72|. 



fluctuations coincides with the correction to the spectrum of 1Z, 






qdS 



6P* 



(72) 



where, using the slow- roll approximation e ^C 1, the spectrum V-r is the one for quasi 
de Sitter space (qdS). This should be taken into account when physically interpreting 
the plots in the subsequent Figs. [4f[T0| because it is the fluctuations in 7Z that are 
directly related to the observable temperature fluctuations (this will be further ex- 



plained in the discussion section, cf. also (34)). The correction to the spectrum SV-r, 
vanishes at the initial hypersurface 770 = — H^ and approaches a non-zero value at 
T) = 0. This means that by the end of inflation an imprint of a small black hole on the 
spectrum will remain. 

For a homogeneous background the propagator in (double) momentum space con- 
tains a delta-peak in the momentum p that is associated with the average of the 



positions. This is seen explicitly in the de Sitter case from Eq. (70). In the case of a 
small inhomogeneity we find a power law divergence at k = q + p/2 = 0. The behavior 
of the spectrum close to this singularity is shown in Fig. [5] 

This observation suggests that an expansion in powers of p ■ q and p, which is the 
basis of the analysis in Ref. [T3] is inappropriate for the complete analysis of small 
black holes in inflation, and, in what follows, we shall not make use of this expansion. 
Instead, we shall analyze the mixed space spectrum without making an expansion in 
powers of p. 

From Eqs. (69) and (67) it follows that i<5G-| = at the initial hypersurface 

(7/ = rjo). In other words, we consider a primordial black hole that was created at a 
time r\ = 770, and we study how it perturbs scalar quantum fluctuations during the 
subsequent inflationary period. 
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Figure 5: In this plot we show the spectrum in (double) momentum space close to the singular 
region for fi = 1/10 and r\ = 0. For this we choose the angle between p and q to be = 0.9ir and 
q = p/2 (black curve, dotted), q = p/2 + Ho (red, dashed) and q = p/2 + 2Ho (blue, solid). We 
observe oscillatory behavior with decaying amplitude. The spectrum diverges for q = —p/2. 



For /j«l also \xp <C k + k' holds and we can expand J_| . |_ in (67) to get 

Jj4 e i(k+k')(r)-r) ) 



■h 



-(vo;p,q, , n,r)) 



l 

kk' 



_JVo_ 
k + k 1 



Akk'(jip) 2 



cos(/ip?7o) 



V 



1 i(k + k')rj 
kk' + kk' 
sm(fipr] ) 



(73) 



fJ-PVo 



1 



»/o 



kk'(k + k') 2 {k + k') 4 (k + k') 3 
///,, T] 2 | irjo \ sin(fiprio) 



cos(npr]o) 



kk'(k + k') (k + k') 2 {k + k') 3 



P-PVo 



where we kept sin and cos unexpanded. From this expression we see that there is a 
contribution that remains finite at future infinity, i.e. at the end of inflation, meaning 
that this correction to the spectrum is propagated through the radiation and matter 
dominated epochs of the Universe. Taking the limit 77 — ^ one finds 



ImJ + _ _+(%;#, Q,V,il) = 

' r] cos((k + k')r]o) sin((fc + k')rj ) 



(74) 



HI 



Ak 2 k' 2 (^p) 2 

+ (W) 2 
1 



k + k' 



cos(fiprj Q ) - 



sin((ipr] ) 



r] cos{{k+k')T] ) 

—3 cos(/ip77 ) 



(k + k' 
2 \ sin((fc + k') m ) 



kk' j \ P-PVo 

1 1 \ ?7ocos((fc + fc , )^o) sm(npT] ) 

kk' (k+k') 2 J (k + k') fipr] 

r]l sin((fc + k')r}o) sin(npr/ ) 



kk' (k + k') 2 J (k + k') 2 
0( V ). 



cos(nprjo) + 



(k + k') 2 



wno 



Note that ImJ^ t ^ (77; p, g*, 77, 77) = 0(rf). Therefore, this contribution to the propa- 
gator is subdominant at late times and becomes completely negligible by the end of 
inflation. Since we are primarily interested in late time cosmology, from now on we 
shall not consider these terms. However, we should keep in mind that these terms be- 
come increasingly important at early times when 77 approaches 770, since they guarantee 
that iSG a i, — > when 77 — > rj a . 
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6.2 Mixed space representation 

Rather than working with the momentum p, we shall mainly consider the mixed space 



propagator by means of a Fourier transformation of (69) with respect to p. The 

mixed space propagator iSG^ {y,Q,Vi r f) ^ s a function of the relative momentum <f, 

the average position y — (x + x')/2 and the times r/ and rf. This makes the physical 
interpretation easier, because the corresponding equal time statistical propagator 

SF(y, q, r), rf) = Re i8G+- {y, g, 77, 77) (75) 

is closely related to the Boltzmann distribution function Sf (6f ~ qSF), cf. [40], and 
hence allows for a simple statistical interpretation as the phase space density of parti- 
cles with momentum (fat position y at time 77. Likewise, the corresponding spectrum 
5V = [q 3 /(2ir 2 )]5F can be given an analogous simple statistical interpretation. The 
dependence on p as well as on y in the mixed space propagator signals a breakdown 
of translational invariance. Furthermore, the mixed space representation is advanta- 
geous also because the power law divergence at p = — 2q, that plagues the double 
momentum space propagator, will become a mild logarithmic divergence in the mixed 
space propagator, whose origin is the IR divergence of the BD vacuum. From the 
above discussion we know that this divergence is regulated when the IR of the de 



Sitter state is made finite. Just as we have regulated the de Sitter propagator (54), 
we shall regulate this divergence by placing the Universe in a comoving box of size 
L = 2n/k . The mixed space propagator can then be written as 

d 3 p 
32tt/j 3 f d 3 p 



iSG + _ (y,q,T],V,k ) = / -—riSG + ^{p,q,r],ri)e lp ' v (76) 

Jk>k ( 27 



I 7 ^3fc /2 (3cos 2 ^-l)ImJ + _,_ + (r 7 o;p,9,ry,r ? )e^. 

Jk>kn \ Z7T ) 



ml J k > k0 (2tt) 

An inspection of this integral shows that it is indeed logarithmically divergent in the 
IR when p = —2q, where k = 0. The divergence is regulated by imposing k > k^. 
The terms in the integrand that cause this logarithmic divergence go as 1/fc 3 for 

small k in ImJ^ , |_ (the integral is convergent at p — 2 q where k' = 0). Let us 

split i5G-\ into an IR finite part and an IR divergent part. For this we choose 

the z-axis to point in the direction of q and keep p general, i.e. q — q(0,0, 1) and 
p = p (cos0sin#, sin^sin^, cos#), and introduce x = cos9 and w — p/(2q). The 
momenta k and k! are then simply 





k = qyl + 2wx + w 2 


ermore, 






3cos 2 0-l = 



k' = q\Jl - 2wx + w 2 . (77) 

3{x 2 - 1) 



1 — 2wx - 



(78) 



The IR limit, k — 0, is then given by the point w — 1 and x = —1 in the (w — x)- 
plane. The other momentum takes the value k! — p = 2q there. So, we can split the 
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propagator as follows 



32tt/j 3 f d 3 p 



(79) 



/r r d 6 f, 

TJW) 



3H t 

32tt/i 3 
~3Hf 



fc' 2 (3cos 2 9 - l)Im J + _._ + ( Vo ;p, q, 77, n)f&* - F d[v ( V , m ,p, q, y) 

F dlv (v,vo,p,q,y), 



Jk>k 



fco ( 27r ) 



with 

F dlv (v,vo,p,Q,y) 



H$e- 2i ** J sin(2 g??0 ) ( ain(2 m ) 



(80) 



M 



sin(2/jg77 ) sin(2g?7 ) 
?7o cos(2<7r7 ) — cos(2(iqr)Q) 



2fiqrj 



2q 



Restricting the momentum k to be above some IR cut-off fco translates to 

■„2 



x> 



"(ShW-^) 



for 1 — ko/q < w < 1 + fco/q. It follows that 

dwdx 



k > ko (1 + 2wx + w 2 ) 3 / 2 



= log 



+ 2 + O(fc ). 



(81) 



(82) 



We can neglect terms that vanish for ko — > 0. As a result, we obtain an explicit 
expression for the cut-off dependent part of the propagator 



32tt/j 3 

3Hq Jk>k 



" P rpdiv/ - - -A 

— 3-F {v,vo,p,q,y) 



o(2-) 
8/iH§e- 2 ^y j shx(2qr]o) 



3nq 2 



2q 



cos(2fiqrj ) - 



sm(2 fj.qr] ) 
2fiqi]„ 



sm(2/iqrio) sm(2qr) ) s 

?/o cos(2qr) ) — cos(2fiqriQ J 



ZmVo 



2q 






(83) 



The power spectrum is defined by (cf. Eq. (52 1 ) 



■a2|o(.r.//).X:P'. v ||0) = / ^.V^y^f^l 



(84) 



with r = ||x — x'||. We have already derived the power spectrum for the de Sitter 
background, equation ( 53 ) . Hence, we will only be interested in the correction induced 
by the black hole, 



8P<t,(q, V, V) = j~2 Re iSG +- (& <?, V, V) ■ 
Again, we can make a split into IR finite and IR divergent part, 
SPM ft V) = SV^ n (q, y, r,) + SVf v (q, y, r,, k ), 



(85) 
(86) 
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with 

-)div 



xr >iivf-> - i. > 4fiqH cos(2q-y) \sm(2q/H ) f sin(2fj,q/H ) 



3tt 3 



2q/H 




fn iv \ am ( 2 w/ H o) sm(2qH ) 

cos(2q/H )— — —j- cos (2fiq/H ) 

2(iq H 2q H 



This expression will be added to the numerical results for the finite part to determine 
the total correction. The finite part is given in terms of an integral in Appendix C, 



Eq. (1161. 



7 Numerical results 



7.1 The anisotropic case: different angles 

For y ^ we can study the dependence of the spectrum on the angle between q and 
y. The numerical results are shown in Fig. [6] All plots show the correction to the 
spectrum SVd, normalized by the scale invariant de Sitter spectrum, Eq. ( 53 ) , 



V, 



<is 



))=0 



m 

4tt 2 



(87) 



Although a good fitting curve has not been found for general angles, it is worth 
noting that besides the high frequency oscillations a modulation with a much lower 
frequency is present which is determined by the mass parameter \i. Therefore, the 
mass of the black hole can, in principle, be inferred directly from the corrections to 
the scale invariant de Sitter spectrum at large q, i.e. at small scales. However, the 
modes of this magnitude that are stretched into today's observable scales are due to 
black holes that were created when the Hubble radius was much smaller than today's 
Hubble radius. The probability to find such black holes inside our Hubble volume was 
estimated in section [31 We will come back to this point in the discussion (see also 
Fig. [111). 



The anisotropic case, y ^ 0, with the vectors yand q being perpendicular (6 = n/2) 
turns out to be identical within numerical precision to the isotropic case, y = 0. 
From the point of view of observations it is hence impossible to associate data to the 
one case or the other. On the level of the integral expressions ( |116| Appendix C) 
for the two cases one can argue that, because of the pole at (w,x) = (1,-1), the 
Bessel function Jo(2qyw\/l — x 2 ) that is present in the anisotropic case is effectively 
evaluated to unity, yielding the same result as the isotropic case. We demonstrate this 
in Fig. [7| Small values of q (large scales) are relevant for black holes that formed at a 
later stage of inflation although their formation probability is exponentially suppressed 
(see discussion section). We present numerical results for the spectrum in the region 
of small qy and q in Fig. p] Furthermore, for q/Ho <1 we can expand 6V$(q, y) and 
find that 



SV^>(q,y,ko) = - ,„ F(qy, 



5„3 



9ir 3 H, 



1) 



16fi b q 
9tt 3 # 



log 



2 I cos(2qycost 
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Figure 6: The correction to the spectrum is plotted as a function of the momentum q for fixed 
qy = 1, fi = f/10 and cut-off ko/Ho = 1/fO but different angles between q and y. Upper left panel: 
8 = 0, upper right: 8 = n/8, lower left: 8 = 7r/4 and lower right: 8 = 7r/4. The case 8 = it/2 is 
treated separately (Fig. [7l. This is the result of a numerical integration with a discretization of eight 
points per unit. For general angles we could not find a good fitting curve. Nevertheless, we observe 
cos(2q/Ho) oscillations modulated with cos(2fiq/ Ho). 
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Figure 7: It is shown that the anisotropic case (blue dots, dark) with vectors q and y being 
perpendicular, 8 = tt/2, coincides with the isotropic case, y = (orange dots, light, mostly covered 
by blue dots). This is illustrated for mass parameter /i = 1/10, cut-off ko/Ho = 1/10 and qy = 1. 
The black fitting curve is given explicitly in Appendix D. 
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Figure 8 : This plot shows the numerical result for the spectrum for small values of the two variables 
q and qy for fi = 1/10 and <(q,y) = 0. A discretization of eight points per unit was used. The 
momentum cut-off is kg/Ho = 1/10 for the left panel and hy/Hy = 1 for the right one. There is a 
significant cut-off dependence. 



with 9 = <(q,y). This means that the cut-off dependent and cut-off independent part 
of the correction to the spectrum factorize individually for low momenta. The integral 
expression of the functions F is given in Appendix C, Eq. ( 118 ). In Fig. [9]we plot the 
corrected spectrum for different angles and small q/H Q . We find that the function F 
is very well approximated by 



F(qy,e) 



-4.8(tog( W ) 

-2.0(\og(qy) 



1.2)cos(2qycos(9)) 
1.3) cos (2qy cos(0)) 



for qy < 0.5, 
for qy > 2.0. 



(89) 



for any angle 9. The fit is not very good in the intermediate region 0.5 < qy < 
2.0. Thus, we have found an analytic expression that describes the correction to the 
spectrum very well for large scales. Moreover, the mass parameter fj, for the black 
hole can be found from the amplitude of the spectrum and the IR cut-off k Q from the 
enveloping curve. 



7.2 The isotropic case: cut-off and ji dependence 

We pointed out in the previous section that the isotropic case virtually coincides with 
the case that the vectors q and y are perpendicular but non-zero. Yet, we can study 
the dependence of the spectrum on fi and the cut-off fcp . The results are summarized 



in Fig. 10 with general fitting curves (c/. Appendix D) that approximate the numerical 
data very well, in particular for small /x. Again, the mass parameter determines the 
modulation frequency. 



8 Discussion and Conclusion 

In this work we have derived the correction to the scale invariant power spectrum of a 
scalar field on de Sitter space from a small primordial black hole to lowest orders in its 
mass parameter /j, = (GMH/2) 1 ^ 3 . To this end, in sectional we have first analyzed 
the probability of black hole formation in the pre-inflationary Universe. In order 
to maximize the formation probability, we have assumed that the pre-inflationary 
Universe is dominated by heavy non-relativistic particles, in which case sub-Hubble 
perturbations grow. We have found that there is a range for the particle mass m < m p 
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Figure 9: The correction to the spectrum as a function of qy is plotted for different angles 9 between 
q and y in the limit q/ Hq — > 0. Upper panels: 9 = 0, middle panels: 9 = 7r/4 and lower panels: 
9 = 7r/2. For all plots we chose q/Ho = 1/10 and fi = 1/10. The upper and middle panels have a cut- 
off kg/Ho = 1/10. The lower panels show also the cut-off dependence: green (upper) dots correspond 
to k /H = 1/100, red (middle) dots to k /H = 1/10 and blue (lower) dots to k /H = 1. The 
fitting curves are given by Eq. l|89| 



and the Hubble rate H* < m p during that period for which the expected number of 
sub-Hubble black holes per Hubble volume can be larger than 0(1), as can be seen in 
Fig.0 

To determine the correction to the spectrum in section [5] we have derived an ana- 
lytic expression for the momentum space propagator of the massless, minimally cou- 
pled, scalar field on the Schwarzschild-de Sitter background in the Schwinger-Keldysh 
formalism. We observe that the propagator diverges in the infrared, demonstrating 
that, in contrast to a recent proposal |13) . an expansion in the momenta is inappro- 
priate, although the breaking of homogeneity is weak. This divergence can be traced 
to the well known infrared divergence of the massless scalar propagator on de Sitter. 
We have used a simple regularization, which consists of placing the Universe in a 
large, but finite, comoving box. In section|4]we have outlined the procedure which has 
allowed to determine the impact of an inflationary black holes on the CMB and struc- 
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Figure 10: We present here the numerical data (red dots) for the correction to the spectrum 
normalized by the de Sitter spectrum in the isotropic case for cut-off ko/Ho = 1 (left panels) and 
ko/Ho = 1/10 (right panels). Upper panels: fj, = 1/10, middle: fi = 1/20 and lower: fj, = 1/50. The 
analytic expression for the fitting function (blue curve) can be found in Appendix D. 



ture formation. We have demonstrated that, to leading order in the dimensionless 
black hole mass parameter /i = (GMH/2) 1 / 3 , the Sasaki-Mukhanov field (curvature 
perturbation) remains the correct gauge invariant, dynamical scalar perturbation. By 
working in the zero curvature gauge, we have then shown how to connect the spectrum 
in the scalar field fluctuation to the spectrum of the comoving curvature perturbation. 
Finally, the knowledge of the appropriate transfer functions has allowed to relate the 
inflationary curvature perturbation to the CMB temperature fluctuations and to the 
large scale structure of the Universe. That analysis is yet to be done in detail. Tem- 
perature fluctuations are obtained from fluctuations in the Sasaki-Mukhanov variable 
by the Sachs- Wolfe effect which yields on large scales 



5T(x,t)5T{x' ,t) 



T 2 



= -(n\n(x,t)TZ(x , ,t)\n) 



(90) 
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The power spectrum V-r of scalar cosmological perturbations was studied in the mixed 
representation in the sections |6.2| and [7] and can thus be related to the temperature 



fluctuations on large scales by performing a Fourier transformation of ( 90 1 with respect 



to r = x — x! . On small scales the more complicated relation (34) has to be used, 
whereby the appropriate transfer functions have to be determined numerically. Other 
transfer functions have to be used if one wants to study the effect on large scale 
structure from perturbations of 1Z. 

The scalar field propagator in the mixed representation is closely related to the 
Wigner function, and hence admits a probabilistic interpretation characterizing the 
Boltzmann distribution function. In section [7] we devoted quite some effort to analyze 
the mixed space spectrum which is a function of not only the relative (comoving) 
momentum q = ||<f||, but also of the comoving black hole distance from us, y = ||y||, 
of the angle <(q,y), and finally of the lowest infrared (cut-off) momentum ko that 
can be excited. Our results are mostly analyzed as the black hole contribution to 
the spectrum relative to the scalar contribution in de Sitter space. Since the observed 
spectrum is highly isotropic, and seemingly homogeneous, we were primarily interested 
in the case when the perturbation induced by a black hole is small, which led us to 



consider the limit /id, or, more precisely, 0.027 < jU -C 1, cf. Eq. (12 1. The lower 
limit comes from the requirement that, before it evaporates, the black hole must last 
at least several e-folds during inflation. The effect of the black hole evaporation during 
inflation is illustrated in Fig. [2j 



The spectrum is first analyzed in section 7.1 for the general anisotropic case, y ^ 0, 
with y the displacement vector of the black hole with respect to us for different angles 
between y and the momentum vector q that is conjugate to the relative distance of two 
points. Then we considered the large scale region (q/Ho and qy small) for the special 
case that q and y are parallel. Furthermore, by making an expansion for q -C H of 
the integral expression for the spectrum we showed that it takes a relatively simple 



form, Eq. (88). We presented an explicit expression for fitting functions for the qy 
dependent part which was determined numerically. In the isotropic case, y = 0, the 
dependence on the parameter fj, is shown. Unlike the spectrum of scalar homogeneous 
perturbations in inflation, which is a function of the momentum magnitude g, and 
depends on two parameters, H and e = —H/H 2 at the Hubble crossing, the scalar 
spectrum in Schwarzschild-de Sitter space depends on q, the distance to the black hole 
y and the angle between q and y. It can serve as a six dimensional template, whereby 
the template parameters are the comoving black hole position y, its mass parameter 
/x, and H and e = —H/H 2 at the first Hubble crossing during inflation. As a summary 
of the numerical results we can conclude that (i) the spectrum as a function of q is 
modulated with a lower frequency which is characterized by the mass parameter /i of 
the black hole; (ii) the enveloping amplitude of the spectrum scales logarithmically 
with the IR cut-off k , and (iii) the isotropic case cannot be distinguished from a 
particular configuration of the anisotropic case (where q-Ly). 

We should point out, however, that a comprehensive analysis of the power spec- 
trum of scalar cosmological perturbations induced by small inflationary black holes 
should take into account also the degrees of freedom of the graviton and start with 
the action for the scalar field and the graviton on the unperturbed background. In 
particular, it is important to study the possible mixing of the scalar, vector and tensor 



modes induced by the black hole. Provided Eq. (32) is an accurate expression for the 



curvature perturbation in terms of the scalar field perturbation in the inhomogencous 



case when a small black hole is present, then based on Eq. (32) one can calculate the 
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spectrum of comoving curvature perturbation induced by an inflationary black hole. 
This is the quantity of interest since it sources the CMB temperature fluctuations 
and the large scale structure. Of course, a more realistic inflationary background is a 
Schwarzschild black hole in a quasi-de Sitter space (SqdS), in which the Hubble rate 
and the deceleration parameter are both slowly varying functions of time. In order 
to make sure that the results presented here hold also for that case, one would have 
to make a complete analysis of cosmological perturbations on a SqdS space, which is 
beyond the scope of the present work. 

Corrections to the power spectrum from primordial black holes are potentially 
relevant for observations. To assess this possibility more carefully, we have to establish 
a relationship between the probability of pre-inflationary black hole formation from 
section [3] and the expected number of black holes that leave an imprint on today's 
CMB sky. A first step in doing that is to trace physical wavelengths, A p h = a/q, back 
in time and check which modes correspond to today's observable scales. In Fig. [TT] it 
is shown how physical scales are stretched during the evolution of the Universe. The 
time axis is given in terms of the scale factor log(a). We are interested in the modes 
which cross for the first time during inflation the Hubble scale at £i x when q = ai x i?i x , 
evolve on super-Hubble scales until some time in matter (or radiation) era, when they 
cross the Hubble scale for the second time at q = a2 X i?2x- At the end of inflation 
the ratio of the physical wave length to the Hubble scale will be, (a/k)/H~ 1 ex a 1_c , 
where we assumed that H ex a~ e and e <C 1 and constant. Since during radiation and 
matter era, H~ x ex a 2 and H^ 1 ex a 3 / 2 , the following relation holds, 

1/2 
«end ^ _ a cq v ( «2x ^ ,„-. n 



&lx / a cnd \ a 



where a eq is the scale factor at radiation-matter equality. Neglecting e <C 1 and taking 



account of the definition, N\ K = log(a on d/ai x ), Eq. (91 1 yields 



A^logfM+ilogf^Y (92) 

\a on d/ 2 V a oq/ 

Now, recalling that a a /a cq ~ 3200, and log(a cq /a cn d) = (1/2) log(-ff e nd/#cq), H Q = 
1.5 x 10~ 42 GeV, H cq ~ 3 x lO" 37 GeV, Eq. ([H| can be recast as 



2 HlO 13 GeV/ 2 
where Z2 X = (ao/a2 X ) — 1 is the redshift at the second Hubble crossing and oq is the 



N ^^ 6l + l h <W^v)-l l ^ 1 + z ^ (93) 



scale factor today, such that the last term in ( 93 1 drops out for modes of the Hubble 
length today \^j 

The black holes that we have so far considered formed during a pre-inflationary era 



(N > 60). From Fig. 11 we see that if a black hole forms between N = 70 and N = 60 
and is within our past light-cone then the scales q/Ho ~ 1 — 100 will correspond 
to observable scales today. Moreover, their amplitude on the CMB will be of the 
order of 10% of the signal obtained for a homogeneous background, as can be seen 
from our numerical results in Figs. [6f |10| Such a strong signal has not been observed, 

12 The absence of the graviton signal in the CMB limits H cnd < 3 X 10 13 GeV (V^V < 1.1 X 
10 16 GeV), and hence ./Vi x ~ 61.6 is an upper limit for the modes with Z2x — 0. In Fig. |ll| jVi x = 60 
e-folds corresponds to about H cnd ~ 10 12 GeV {V^{* ~2x 10 15 GeV). 
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log(scale) 



toservable 
scales 



f = 2GM 




N = 60 



N = 



today 



Figure 11: The stretching of modes during the evolution of the Universe is displayed. The time 
axis is log(a/a cnt i) = — N with N the number of e-foldings before the end of inflation. During matter 
domination the Hubble rate is given by H(£) oc a -3 ' 2 , during the radiation era it is H(t) ~ a -2 . A 
constant Hubble rate Ho during inflation (slow-roll approximation) is assumed. The upper diagonal 
line (red) shows the stretching of the mode that re-enters the Hubble volume today, thereby defining 
the largest observable scale. That mode crossed the horizon approximately No ~ 60 e-foldings before 
the end of inflation (sometimes a value of Nq = 50 is assumed). The smallest observable scales depend 
on the maximum resolution of the apparatus used for the measurement. Today's best resolution is 
/ ~ 7r/#(rad) ~ 10 2 . The lower diagonal line (red) marks this mode, it exited the Hubble horizon at 
N = 50. 



implying that such black holes had not formed. If inflation lasted much longer than 
60 e-foldings, it would be useful to estimate the probability for black hole formation 
during inflation. In order to estimate the formation probability in the early stages of 



inflation, we first note that on sub-Hubble scales Eq. ( 14 ) for density perturbations 
can be used in inflation Fj with the equation of state parameter w ~ for matter 
and the Hubble rate H ~ Hoa~ e . In slow-roll inflation we can take the Hubble rate 
to be constant to estimate the growth of the perturbations. Since po(t) ~ l/a 3 ( 1+w \ 
the last term in ( 14 1 becomes rapidly negligible and we find that the evolution of the 



perturbation is given by 



6{k,t) = 6*(k) + 6.(k)e 



-2fHdt 



(94) 



which means that the leading perturbation tends to a constant. Using the relation 
between the density perturbation and the fluctuation in the number of particles in the 



13 cf. also Lyth and Liddle 'The primordial density perturbation'. 
(2009) 



Cambridge University Press 
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comoving ball with physical radius R, Eq. ( |T8| , it follows that during inflation (t > ti) 

5N(R,t) SN(R,h) 



(N(R,t)) (N(R,tj)) 



(95) 



where tj denotes the beginning of inflation. Therefore, also the variance in the fluc- 
tuations of N(R, t) is constant for t > t], 

a(R,t) = (SN 2 (R,t)) = J^0r (N(R,t)) 2 = (6N 2 (R, tl )) (96) 

because the average number of particles in a comoving ball does not change during 
inflation p^ S ince the universe is dominated by the potential energy of the inflaton 
field, Eq. ( |21[ ) is not valid during inflation. Taking this into account, the Friedmann 
equation is 



Stt 



H S? W)+VW) ' (97) 

and we find that the critical fluctuation that is necessary for black hole formation is 
given by SN CI (R,t) = rrvLR/(2m) — (N(R,t)} sw miR/{2m). Considering fluctuations 
in a volume of radius R, we thus obtain the result 



5N CT (R,t) m 2 P R 



(98) 



y/2a{R,t) 2m v / 2(SN 2 {R,t I )} ' 
In the case when R is the comoving radius, R = Rze Nl ~ N ^\ the expression in (98) 



grows exponentially in physical time. On the other hand, when R is a constant 
physical radius, yields 8N CT /y/2a oc e?( Nl ~ N w\ We conclude from Eq. (25) that the 
probability for black hole formation is rapidly (faster than exponentially) suppressed 
during inflation. That means that if inflation lasts much longer than about 70 e- 
foldings the probability of observing any inflationary black holes within our past light- 
cone will be tiny. Furthermore, from Figs. |9| and [TT| we see that the black holes which 
form later during inflation (corresponding to N < 50) are observed as q <C Hq and 
induce unobservably small fluctuations in the CMB. According to the above analysis 
the probability for formation of these black holes is tiny. 

We shall now comment on various ways how this work can be extended. The 
case of a black hole in inflation can be easily generalized to heavy particles and in 
particular monopoles by gluing the exterior Schwarzschild-de Sitter geometry to an 
interior matter solution. The possibility of magnetically charged black holes in the 
presence of magnetic monopoles is discussed in |41j . Also, it would be interesting to 
study a rotating body giving rise to an exterior Kerr-de Sitter geometry. Furthermore, 
there is no reason why a black hole should form in the rest frame of the inflaton 
field. Assuming a black hole forms with an initial relative speed vq, its motion will 
be damped the same way as the motion of a test particle, which obeys a geodesic 
equation, v + Hv sw 0, such that the speed is Hubble damped as v = vq/cl after the 
formation. Nevertheless, in the early stages after formation the speed of the black 
hole can yield significant velocity perturbations which may have observational impact 



14 To see this consider the geodesic equation for particles during the inflationary phase, v + Hvzx 0, 
which shows that their motion (with respect to the fluid rest frame) is exponentially damped, v(t) = 
#(£r)e — ( "i~ "(*)). This means that very quickly particles get frozen and cannot leave the comoving 
volume, leading to a constant particle number per comoving volume. 



M2 



on the structure formation, and an interesting question is whether a relation can 
be established with the recent observation of large scale flows [T5]. It is also worth 
investigating what is the damping of the speed of a spinning massive object. 

When this work was nearing completion, we became aware of Refs. [J5J H3J HI] 
which, just like this work, address the problem of scalar field fluctuations on a Schwarz- 
schild-de Sitter (SdS) background. We note however that there are important differ- 
ences between these papers and ours, and that our analysis of the black hole per- 
turbations on SdS background is much more elaborate [j£| Since the first version 
of the present paper, several articles have appeared which discuss anomalies in the 
CMB from black holes and point-like defects. In a recent proposal by Gurzadyan 
and Penrose [JS] US] the attempt was made to relate the presence of (anomalous) 
concentric circles in the CMB sky to pre-Big Bang activities involving black holes in 
the framework of conformal cyclic cosmology. This has been the subject of recent 
debates [47j [48] . The role of the presence of pre-inflationary particles for anomalies 
in the CMB was investigated in [55J and [53J where the imprint is also expected to be 
given in form of rings. Finally, the effect of pre-inflationary black holes on the CMB 
power spectrum was studied in |49] and a link to the quadrupole anomaly was made. 
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Appendix A 

The coordinate transformations that give rise to the line element in cosmological form, 
Eq. Q, are presented here. 

In its static form the Schwarzschild de Sitter (SdS) solution is given by the line 
element 

ds 2 = -f(r)dl? + f(r)- 1 dr 2 +r 2 dn 2 (99) 

/(f) = 1 - 2GM/f - f 2 /R 2 , 

with R — y 3/A being the Hubble radius and A the cosmological constant. This line 
element reflects the spherical symmetry and the time translation symmetry. 
If we take 

T(t,r)=t-J v 1 ~H r) dr, R{T,r)^f/a{ T ) (100) 

for a(r) = e T / R ° then the metric becomes 

ds 2 = -dr 2 + a 2 {r) [{dR + F{t, R)dr) 2 + R 2 dfl 2 ] (101) 

15 The most important difference is that, in contrast to this work, Cho, Ng and Wang 42 do not 
make any attempt to relate the inflaton fluctuation to the comoving curvature perturbation, which 
can be defined on the Schwarzschild-quasi de Sitter space, and which is the correct generalization of 
homogeneous inflationary spaces. Furthermore, we provide an analytic estimate for the formation 
probability of small pre-inflationary and inflationary black holes. Next, we provide an exact answer 
for the Schwinger-Keldysh propagators in the double momentum space (whereby admittedly treating 
the black hole as a small perturbation on de Sitter space), and we provide a more complete analysis of 
the significance of the results by discussing how the curvature spectrum gets perturbed as a function 
of the position of the black hole. 



:W 



with 



R 



hGMRl 



F ^ R) =R- [ l -VR^(r) 



1 



(102) 



The spatial slices r = const, are flat. In order to make the metric diagonal, we have 
to solve a differential equation for a function R(t, r) 



dR 



-F(t,R), 



which has the solution 



R(r,r) = {2GMR 2 Q ) 1/3 e~ T/Ro smh 2/3 



3 /^3 r + T (r) 
Ru 



(103) 



(104) 



where To(r) is a r-independent integration constant. This constant can be determined 
by the requirement that R(r,r) ~ r for M = 0, i.e. homogeneous cosmology is recov- 
ered. One finds that 



T (r) = i? log 



GMRl 



1/3 



(105) 



and hence, 



. GMRr. 



2\ 2 /3 



The metric then becomes 
dR. 



ds 



-dr 2 + a 2 {r) 



dr 



(r,r)) dr 2 + R 2 (T,r)dn 2 



(106) 



-dT 2 



(2GMR 2 ,) 



2^2/3 



dnh 4 / 3 



3r + T (r) 



^0 



( . o);h 2^3 r + To(r) ]</f , ): , /()J 
2 -Ro 



-^^^-§^'" 



GM_Rg 
2a(-r) 3 r 3 

GMR 2 
2a(-r) 3 r 3 



dr 2 



w 



which is recast in the main text, Eq. (pi), in conformal time rj, dr\ — dt/a{t). Note that 
the metric is singular at 



t (t) 



( GMR 2 \ 1/3 1 



(107) 



and regular for all values r > ro(r). As a check, R^ — Ag^ — 0, and the Kretschmann 
invariant is 



p upvpa _ 

± ^pvpo ± ^ 



48GM 2 8A 2 



48GM 2 



r.{) 



«(^ 6 (l-2Al^) 4 



8A 2 



(108) 



with f denoting as previously the Schwarzschild radial coordinate. This proves that 
ro corresponds to the curvature singularity at f = 0. Since df/dr > for r > ro(r), 
we conclude that f —$■ oo when r — > oo without any double-valued regions. 
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Appendix B 



For the derivation of ( 62 ) from ( 60 1 we have to evaluate first 
dr" rl 



Ii{k,k',r)") 



h(k,k',r)") = 



dcos9" 

i ./D 



y e i(fc'-fc)-£'' 



4;, I -Ci(p) 



(109) 



' ; dcose" / #" 



4irk 12 



\k ■ x ) j(k'-k).x" 



r ll2 



smp 



- Ci(p) 



3 cos 2 6»-l 



cos p sin p 



P* 



P° 



with the cosine integral function Ci(z) = — J dtcos(t)/t, 9 — <(k' — k,k') and 
p = -p\\k-k'\\r]". Then, it follows that 



I(k,k',r)") = 



oc 



dr 



n r i 



,■" 



2tt / 

dcos9" / d<b" \ -k' 2 + 3 



fj,r)" ■ J -\ JO 

-k' 2 h{%,k', rf') + 3I 2 (k,k',r)") 



(k' ■ x") 2 

r ll2 



e i{k'-k)-3" ( U0 ) 



Auk' 2 (3 cos 2 6-l\ 



cos(/i||/c - k'\\r)") sm(p\\k - k'\\r)") 
/_t 2 ||fc-fc'||V 2 ^ 3 ||fc-fc'|| 3 ?7" 3 



To establish (64) and d65fo we used the following relations for the step function: 



9(77 - v")®{v" - rf) = @(v - rf) [Q(rf' - rf) - @{r)" - v)} (HI) 

Q{r) - tf')Q(rf - rf') = ©(»» - rf)Q(rf - rf') + ©(»/ - rj)&(v - »/') 

6(7?" - 77)9(77" - r/) = 8(77 - 7/)9(7/' - rj) + 6(77' - 77)9(77" - 77') 
9(77" - 77)9(7?' - 77") = 9(77' - 7;) [9(7?" - 77) - 9(7/' - 7/)] . 



The final form for J^ _| and J_| [_, Eqs. (67) and (68), can be obtained by solving 



the following integral, for real parameters A, B, a and a ^ 1 and p = —prj"p with 
A + aB = 0, 



dp (p 3 + iAp 2 + Bp) 



cos p sin p 
p 2 p 3 



a iap 



(a 2 - l) 2 



a(A — ip) 
a 2 -I 



cosp - 



i(a — A) — p 



2ia 



B 



y 2 - 1 (a 2 - l) 2 



Slll/9 



(112) 



^iap 



Appendix C 



We present here the general integral expression for the finite part of the spectrum. 
For this we write 



9 = 9(0,0,1) 

p = p(cos <fi sin 9, sin <p sin 9, cos 9) 

y = y (cos 4> y sin 9 y , sin <f> y sin 9 y , cos 9 y ) 



(113) 
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Introducing x = cos 9 and w — p/(2q), as in the main text, we find that 



/ d<p cos(p '■ y) ~ 2n Jo(2qyw\' 1 — x 2 sin 9y) cos(2qywx cos 9 y ). 
Jo 

In terms of the variables k and k' , 

K = k/q = v 1 + 2wx + w 2 , k! = k' j q = yl — 2wx + w 2 



the integral form of the spectrum, in Eqs. (|T9j) and (85 86), becomes 

3(a; 2 -l)' 



5V^{q,y,rf) 



2/j,qH 



3^ 3 ,o 
cos((k + K,')q/H Q ) 



00 , f 1 dx . 
dw -^2 



(114) 



(115) 



(116) 



K + k' 

cos((k + K r )q/Ho 



cos(2fiwq/H ) 



sm(2[iwq/HQ 



(2/iW) [ (* + *>)* 

2sm((K + K,')q/H ) 



cos(2[iwq/ H[ 



2fiwq/Ha 

sin(2 fiwq/Ho) 



(q/H )( K + K'Y 



cos(2[iwq/Ho) 



2/iwq/Ho 
(q/H ) sin((« + n')q/H ) sm(2fj,wq/H ) 



(k + k') 5 



2/iwq/H 



x Jo(2qyw\/l — x 2 sin 9 y ) cos(2qywx cos 9 y ) 
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sin((K + n')q/H ) 



2fiqH f°° , f 1 dx 



3(x 2 - 1) 



,'2 



(q/H )n' 

(2/iw) 2 



cos(2/iwq/Ho) 



sm(2fiwq/ Hq) 
2/iwq/Ho 



cos((k + K')q/H ) sin(2^wq/H ) sin((« + K')q/H Q ) 

7 ^cos(2fj,wq/H a ) 



k'{k + k') 2fiwq/H Q (q/H Q )n'(K + k') 

x J {2qyw\/l — x 2 sin 0„) cos(2qywx cos # a ) 
JsinjVgo) f cos(WgQ) sin(2 M /i?o) 



2g/#o 



2fiq/H 



n 



to /l7 ^ in ( 2 w/ ff o) sm(2q/H ) , . . 
cos(2q/H a ) — — cos(2fiq/H a ) 



2fiq/H Q 
For g/Ho "C 1 one finds 



2«/fl" 



cos(2qy cos 9y) 



16/iV 
9tt 3 H 



16/*V 

9n 3 H 



SP*(q,v,v) = ^^F(qy,0y) + Sr ( lo s ( s ) + 2 ) cos(2 W cos0 s ) (117) 



:-!(> 



with 
F(qy,6 y ) = 



(118) 



-!■ / diu to 



c/.i; 



x k 2 (k + k') 3 



x 2 + 



3(x 2 - 1) 



2Jo(2qyw\/l — x 2 sm.6 y ) cos(2qywx cos 6 y ) 



+ dw 

Jo 



1 dx 

. ^3 



3(ir 2 — 1)\ 2w A Jo(2qyw\/l — x 2 sm6 y ) cos(2qywx cos9 y ) 



..12 



k'(k + k') 



cos(2qycos8 y ) 



Appendix D 

This appendix contains the fitting functions for the numerical data. In the isotropic 
case y = 0, the correction to the spectrum normalized by the scale invariant de Sitter 
spectrum Aw 2 /Hfi has been fitted in Fig. 10 with the following function: 



f(q, n, ko) = iiA{ji) log(q/H ) [(1 - cos(2 m /H )) sin(2q/H - 1/4) 
+ ^B^iq/Ho) cos(2q/H a - tt/4) 
, 16fiq/H \ sm(2q/H - 1/4) 



(119) 



3?r 



2q/H 



cos(2fj,q/Ho) — 



sm(2^q/H ) 
2 m /H 



i n sm(2uq/Ho) s'm(2q/H - 1/4) 
cos(2q/H Q - 1/4) v m "' v u " '—!■ cos(2 W /5 ) 



x (log(3g 2 ) + 2^ 
16nq/H Q j sm(2q/H ) 



3tt 



2<z/5 



2fjLq/H 



cos(2[iq/Ho) 



2q/H 

sm(2nq/H ) 
2fiq/H Q 



io iv N sin ( 2 M9/^o) sin(2g/5 ) ,_ .„ / 
cos(2q/H ) — — cos(2fiq/H ) 



2 M /H 



2q/H a 






We have no analytic expression for the functions A and 5 but we observe that the 
dependence on (i is weak, 

,4(1/50) = 0.76, 5(1/50) = -1.0. (120) 

,4(1/20) = 0.87, 5(1/20) = -0.7, 
,4(1/10) = 0.80, 5(1/10) = -1.3, 

Obviously, setting A = 0.8 and 5 = —1.0 for all /i results in reasonable fits, too. 
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